Superconductivity in Na2;Co02 yH20 : 
Is Spin-Charge Separation Protecting a di+id2 State ? 
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Superconductivity in Nai:Co02-yH20 is likely to be a p or d-wave; however, experiments are 
unable to pinpoint the symmetry. A simple estimate of pair breaking effects from an unavoidable 
'Na+ vacancy disorder' in an ordered Na+ lattice, at an optimal x^pt ~ 0.30 is shown to destroy 
a Fermi liquid based p or d-wave superconductivity. However, a robustness of superconducting 
and normal states, seen in experiments is pointed out and argued to imply presence of a 'quantum 
protectorate', possibly a 'spin-charge decoupling' that protects a di+id2 and not a p-state. A 
calculation of Knight shift and tjt- in the framework of RVB mean field theory and a fit to the data 
of Kobayashi et al. [9] is made. 

PACS Numbers: 74.20.Mn, 74.20.Rp, 76.60.Cq, 76.60.-k 



The recent discovery of superconductivity (SC) in 
Na2;Co02-yH20 [1] has generated renewed interest, due 
to the role of strong electronic correlations, in the cobalt 
oxide family. This compound has layered crystal struc- 
ture where Co atoms form a two dimensional triangular 
lattice separated by Na+ ions and H2O molecules. A 
number of experimental results [2-14] in these family of 
material are quite anomalous and remains unexplained. 

A problem of current interest is the symmetry of the 
superconducting order parameter. In known examples, 
where electron correlation provides a superconductivity 
mechanism, the order parameter symmetry is unconven- 
tional, either d or p-wave. Nuclear spin resonance experi- 
ments, which can identify order parameter symmetry, has 
so far not been decisive on this issue for Na2,Co02-yH20. 
On the theoretical side, a single band t-J model intro- 
duced by one of us [18] allows, at small dopings, a spin- 
charge decoupled normal state and a resonating valence 
bond (RVB) superconductivity with di-|-id2 symmetry; 
and a triplet pairing at high dopings, arising from a ring 
exchange processes and a consequent ferromagnetism in- 
duced by dopant dynamics. Subsequent papers [19-21] 
have focussed on RVB theory of d-wave pairing and oth- 
ers [22-24] on triplet pairing invoking certain ferromag- 
netic tendency. 

The one band t-J model physics introduced in refer- 
ence [18] implies a distinct possibility of 'spin-charge de- 
coupling' for a wide range of doping. An important con- 
sequence, as we have learned from our study of cuprates, 
is a remarkable possibility of a 'quantum protectorate' 
[26,27]. That is, an ability to protect superconduct- 
ing and non-Fermi liquid normal states from phonons 
and disorder. This issue needs to be investigated in 
Naa;Co02-yH20. However, a quantum protectorate as- 
pect is visible in some experiments; i) replacement of 
H2O by D2O, that alters the phonon frequency of ice 
layers substantially, in the experiments of Jin et al. [5] 
has no appreciable effect either on normal state resistiv- 



ity or the superconducting Tc ii) survival of supercon- 
ductivity in a narrow dome {j < x < ^) [7], overcoming 
a possible strong competition from charge order [25], a 
charge localization tendency, iii) unavoidable disorder in 
the Na-H20 layers [13,14], another charge localization 
tendency. A large ab-plane residual resistivity, pab ^ 300 
lifl cm is observed in a first single crystal (Tc ^ 5 K) 
measurement [5]; even if a fraction of it were intrinsic, a 
Fermi liquid based d or p-wave would not have survived. 
For example, in Sr2Ru04, a well established p-wave su- 
perconductor [28] even a small disorder corresponding to 
a residual resistance of pab ~ 1-5 /xfi cm is capable of 
suppressing a superconducting Tc (« 1.5 K) to zero. 

After summarizing known phenomenological evidence 
for p or d-wave superconductivity, we argue for the pres- 
ence of an unavoidable Na'^ vacancy disorder in an or- 
dered Na"*" plane at optimal doping x — Xopt ~ 0.30. We 
quantify the disorder and show that it is strong enough 
to destroy a Fermi liquid based p or d-wave supercon- 
ductivity. We then discuss how spin-charge decoupling is 
capable of protecting a d-wave and not a p-wave from the 
singular effects of disorder known from Anderson's theo- 
rem of dirty superconductors. A simple RVB mean field 
theory results for Knight shift and ^ and comparison 
with experimental results of Kobayashi et al. [9] is 

TABLE I. Summary of the Knight shift and the nuclear 
spin relaxation rate for Naa;Co02-yH20: 
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made. 

Various measurements have looked for unconventional 
superconductivity in Na2;Co02-yH20: uniform spin sus- 
ceptibility Xs, NMR/NQR, specific heat and resistivity. 
A summary of nuclear spin resonance results [9-12] are 
shown in table 1. Some of the marked differences in 
NMR/NQR results are likely to arise from sample qual- 
ity. Uniform spin susceptibility measurements from all 
the groups [3-6,9,11,12] has the form [3]: Xs ~ XoiT) + 
[T+B) ■ Here Xa{T) is a nearly temperature independent 
but enhanced Pauli component (enhanced by a factor of 
5-7 compared to free electron value); this component also 
exhibits a 'spin gap' like reduction by about 5-10% from 
250 K to Tc. The Cm-ic- Weiss component has a large 
variation between groups: single crystal measurements 
by two groups [5,6] show a small value of Curie constant 
C ~ 0.05 X 10~^ emu K/mole compared to ~ 1.0 x 10~^ 
emu K/mole given by other groups from powder samples 
[3,4,12]. It is likely, as suggested in reference [3,9] that 
the Curie- Weiss component may not be intrinsic at opti- 
mal doping in NaxCo02-yH20. 

Nuclear spin relaxation scales with Xs [12], for a 
range of temperatures above Tc- However, a Curie- Weiss 
enhancement in ^^^^ below about 100 K is not seen 
by all groups. Notably, Kobayashi et al. [9] observe a 
« 12.5 sec-i K-i in the range Tc < T < 30 K 
and no significant enhancement. This result also ques- 
tions the intrinsic character of Curie- Weiss enhancement 

m jTj:. 

Specific heat coefiicient 7 shows, in all experiments 
[3-6,15-17], an enhancement over the free electron value 
suggesting a mass enhancement by a factor about 7. One 
of the single crystal measurements of 7 and Xs [6] esti- 
mate a Wilson ratio R= ^7777^^7- « 1.53. Accordingly 
they interpret the system as a heavy metal. 

Putting measurements for T > Tc together we may say 
that there is no particular preference between a p or a 
d-wave, even though many experimental groups suspect 
a p-wave. 

The situation below Tc is unclear, with conflicting re- 
sults for coherence peak and additional Yoshida-function 
like Knight shift. However, it is fair to point out that 
Kobayashi et a/.'s result is intriguing and demands a fur- 
ther discussion; i) in addition to absence of Curie- Weiss 
enhancement above Tc, alluded to earlier, there is a clear 
coherence peak below Tc in ^ and ii) one sees an addi- 
tional Knight shift and a Yoshida function type behavior 
below Tc- Is it possible that Curie- Weiss contributions 
are extrinsic and Kobayashi et aL's sample has less of it 
and therefore has managed bring out a coherence peak 
and non-trivial Knight shift ? It is important to investi- 
gate this issue further. 

At the present stage it is difficult to conclude about 
possible nodal character of the order parameter from the 
power law behavior of ^ [11] and specific heat [16,17] 



below Tc, until we exclude extrinsic effects. 

Putting measurements for T < Tc together, we may 
say that perhaps there is a support for a gaplful di-|-id2 
or pi+ip2 pairing, imlcss the coherence peak and addi- 
tional Knight shift measurements of Kobayashi et al. [9] 
are proved wrong. 

Triplet superconductors such as, Sr2Ru04, UGe2 etc. 
are known to be ultrasensitive to non-magnetic (scalar) 
disorder. In what follows, we make an estimate of 
pairbreaking effects from two sources of disorder in 
Naa;Co02-yH20 at optimal doping, within the context 
of Fermi liquid theory. We find, contrary to the exper- 
imental observation, that existing pair breaking effects 
should suppress superconducting Tc completely. 

In the theory of dirty superconductors [29] , the param- 
eter a = ^f.^j. ^ plays an important role; here ^ is the 
'pair breaking' rate and ^ is the energy spread an cigcn 
state (a Bloch state) , close to the Fermi level, acquires on 
introduction of disorder in the system. When this energy 
spread is nearly equal to fcsTco, i.e., awl superconduc- 
tivity disappears. This also roughly corresponds to the 
electron mean free path I at the Fermi level becoming 
nearly equal to the coherence length ^. 

Sources of pair breaking are i) an irreducible Na"'' va- 
cant disorder at 'optimal doping', ii) low frequency charge 
order fluctuations and iii) dissipation from the low fre- 
quency polar modes of the H2O layers. We first discuss 
the effect of disorder in Na"*" layer. As elaborated in [25] 
it is convenient to define a reference 'stoichiometric' com- 
pound NaiCo02 • IH2O, where Na atoms arc ordered; 
in an ideal ordered structure they fill one of the three sub 
lattices of a triangular lattice of preferred sites and have 
a \/3 X VS order. This is consistent with the large co- 
herence length of 1000 A in the ordered Na+ layer seen 
in the structural study [14]. Experimentally optimal Tc 
« 5 is reached at Xopt = ^(3: + 5) ~ 0.3; that is, a re- 
moval of (5x = (i — Xopi) w ^ of Na atoms from the 
ordered Na"^ layer. Assuming that these vacancies are 
disordered, we get an 'irreducible' vacancy disorder in an 
otherwise ordered lattice of Na"*" ions. A corresponding 
missing 'dipole moments' exists in the H2O layer as well. 

At optimal doping, a fraction 2 x ^ of vacant Na+ 
sites from the sandwitching Na+ layers above and below 
provide a random potential to the delocalized carriers in 
the C0O2 planes. A change in potential at the nearest 
Co site, arising from a removal of one Na+ is « ^^/^ ^ ^ « 

80 meV; here R w 5.3 A and etotai = ^cc + (H2O ~ 30. 
In the above we have taken into account a large dielec- 
tric screening from the hydrogen bonded H2O sheet [25]. 
The strength of the above scattering potential 80 mcV is 
the same order of magnitude as the bandwidth 200 meV 
as estimated in reference [6]. This makes the mean free 
path of an electron in the above tight binding band of 
the order of the average vacancy distance ~ 4 (in 
units of Co-Co distance). 
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The above estimate of the mean free path is less than 
100 A), the coherence length in our superconductor. 
Thus we get the result that a < acriticai ~ 1, making 
a Fermi liquid based p or d-wavc superconductivity un- 
likely. Since the Na+ vacancy disorder is strong enough 
to suppress a d or p-wave superconducting Tc, we do not 
discuss our estimate of additional pair breaking effects 
from charge order fluctuations and low frequency polar 
fluctuations of ice layers. 

The fact that the superconductor has survived severe 
pair breaking effects of disorder discussed above and a ro- 
bustness with respect to replacement of H2O by D2O [5] 
suggests to us a possible existence of a 'quantum protec- 
torate'. If our 2D t-J modeling is correct, a spin-charge 
decoupled state is a likely consequence and it will be- 
come a natural quantum protectorate. Spin-charge de- 
coupling, though suspected, has not been experimentally 
investigated so far in NaxCo02-yH20, a strongly corre- 
lated 2D metal. The issue of quantum protectorate is 
well known in cuprates [26], another strongly correlated 
2D metal, through the robustness of the non- Fermi liquid 
normal state and d^2_y2-wave superconductivity against 
off-plane disorder and phonons. 

Below we give a physical picture of spin-charge sepa- 
rated state as a quantum protectorate and how a singlet 
d-wave and not a triplet p-wave will be protected. First 
we recall how non-magnetic disorder affects pairing in a 
Fermi liquid: scattering of individual k states of elec- 
trons by scalar disorder interferes destructively the co- 
herent scattering of (k, — k) pair over the Fermi surface, 
whenever it is trying to establish a pair amplitude that 
changes sign along the Fermi surface, X]k^(^) ~ ^ 
in the case of d or p-wavc. 

In a t-J model, the singlet stabilising collision between 
two electrons is a superexchange process that happens 
between two singly occupied sites (neutral spins), in our 
case between two adjacent spin-half moments Co^+'s. 
The charge 2e condensation, namely superconductivity in 
RVB theory is viewed as taking place in two steps: i) pair- 
ing between neutral spins in the spin singlet channel and 
then condensation of holon or doublons. The quantum 
protectorate hypothesis essentially states that the effect 
of scalar disorder on spinon/holon/doublon (collective 
'topological' excitations) is in some sense 'averaged out', 
and remain unaffected for small disorder. To this ex- 
tent a local superexchange process that involves two neu- 
tral spins remains unaffected. Similiarly a holon/doublon 
condensation that involves a spinless charge also remain 
unaffected. From this point of view, as long as the pairing 
is in spin-singlet channel it is compatible with the robust 
superexchange process. Thus singlet pairing such as an 
extended-s or d-wave can be protected by spin-charge de- 
coupling. 

Next important question is how a p-wave may not be 
protected by spin-charge decoupling. A possible p-wave 
pairing suggested in reference [18], goes beyond superex- 



change and invokes a ring exchange [30] involving a mini- 
mum of four sites and 5 electrons and a consequent ferro- 
magnetic tendency given by an effective Jeff « J - < 
(when the hopping matrix clement t has a 'wrong' sign). 
This dopant induced dynamics involves charge transport. 
To this extent it is not a triplet pairing between two neu- 
tral spins. It is rather a pairing between two electrons in 
the triplet channel. Thus spin-charge decoupling is un- 
able to protect a triplet state, that involves fragile elec- 
trons rather than the robust neutral spins or doublons. 

Having offered some phenomenological and micro- 
scopic support for di-|-id2 pairing, we have performed 
a simple RVB di+id2 mean field theory calculation for 
and Knight shift and fit to the experimental results 
of Kobayashi etal. [9]. We point out that Bang et al. 
[31] have calculated the effects of impurities on ^ in the 
di-|-id2 phase. We refer the readers to reference [18-21] 
for details of RVB mean fleld theory. After a quantitative 
fit to Kobayashi et aL's data, we get a value of hypcrfinc 
contact coupling to be 39.32 X 10~'*cm~^. Using the same 
hyperflne coupling, we compute the temperature depen- 
dent nuclear spin relaxation rate and make quantitative 
comparison with the existing experimental data. 

Nuclear spin (Iq) relaxes by coupling with the local 
conduction electron spin density (So) in a metal. The in- 
teraction Hamiltonian is given as, iJhyp = ^Iq-So, where, 
A is the effective hyperflne coupling. The spin part of the 
Knight shift is, 

where, Xs((;, w) is the dynamic spin susceptibility. The 
temperature dependent order parameter in a RVB SC 
is computed self consistently and is incorporated in the 
susceptibility to get the temperature dependent Knight 
shift. We flx the chemical potential /Lt=0.075 eV to get 
30% electron doping. We also take J=6.2 meV in or- 
der to get a Tc=4.5 K. Throughout the computation, we 
have chosen teff=0.03 eV and a cut-off energy scale to 
be 500 K. The ratio 2A(T = 0)/kBTc comes out to be 
approximately 3.5. 

In Fig. 1 we show the flt of the computed temperature 
dependent Knight shift to the data of Kobayashi et al. 
[9] from which we extracted the value of the hyperflne 
coupling parameter, A=39.32x 10~^ cm~^. 

In the process of fitting the data we subtracted a fairly 
large temperature independent backgromid {Korb%) of 
2.55% from Ktot%, where Ktot{T) = Korb + K,{T). And 
Ks{T = Te) - Ks{T = 0) is nearly 0.6%. The abso- 
lute value of Ktot% in this case is nearly 3.15%. The 
value of A estimated as above compares well with high 
Tc cuprates [33] whereas Ktot{T) and Korb do not. In 
cuprates, Korb was estimated to be 0.23% and KgiT = 
Tc) — Ks(T = 0) = 0.37% which are much smaller com- 
pared to layered cobaltate. 
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FIG. 1. Fit of the computed Knight shift in a RVB SO with 

di+id2 gap symmetry (soHd Uncs) and a BCS SC with s-wave 
gap (dashed lines) to the data by Kobayashi et al. [9]. 



Fig. 2 shows that the overall shape of the spectrum 
looks similar to the data of Kobayashi et al. [9] and Waki 
et al. [10] but with some differences. We find that the 
value of Kobayashi et al. [9] is about 0.65 

times the computed value in the present formulation. 
The smallness of the coherence peak seen in the exper- 
iments is likely to be an effect of strong correlations in 
addition to disorder [31,32]. 

To conclude, it will be important to perform exper- 
iments to look for signals for a spin-charge decoupled 
state, a quantum protectorate and chiral di-|-id2 super- 
conductivity. 

The authors would like to thank Prof. K. Yoshimura 
and Prof. Y. Kobayashi for providing NMR/NQR data. 
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FIG. 2. Computed nuclear spin relaxation in a RVB SC 
with di-|-id2 gap symmetry (solid lines) (multiplied by 0.65) 
and a BCS SC with s-wave gap (dashed lines) compared to 
the data by Kobayashi et al. [9]. 

The relaxation rate from the golden rule is given by, 
1 27r 

= ^ \''J2i^ik,'l)fk{i-fk+,)S{Ek+g-Ek-fiwo), 

k,q 

where, the coherence factor P in a di+id2 superconduc- 
tor is written as, l^(k. q) = i (i + £Mi±a+J^Ml^ii±ai). Here, 
Wo is the nuclear Zeeman frequency which is extremely 
small on the scale of energy of interest. In computing 
the relaxation time for a RVB superconductor, we con- 
sider the same parameters as we used for computing the 
Knight shift, also the same temperature dependence of 
the SC order parameter. Further, the hyperfine coupling 
which we estimated by fitting our data with the Knight 
shift measurement has been used to compute the relax- 
ation time. 
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